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Abstract: The properties of planar ice crystals settHng horizontally have been investigated using a vertically-pointing Doppler lidar. 
Strong specular reflections were observed from their oriented basal facets, identified by comparison with a second lidar pointing 
4° from zenith. Analysis of 17 months of continuous high-resolution observations reveal that these pristine crystals are frequently 
observed in ice falling from mid-level mixed-phase layer clouds (85% of the time for layers at — 15°C). Detailed analysis of a case 
study indicates that the crystals are nucleated and grow rapidly within the supercooled layer, then fall out, forming well-defined 
layers of specular reflection. From the lidar alone the fraction of oriented crystals cannot be quantified, but polarimetric radar 
measurements confirmed that a substantial fraction of the crystal population was well oriented. As the crystals fall into subsaturated 
air, specular reflection is observed to switch off as the crystal faces become rounded and lose their faceted structure. 
Specular reflection in ice falling from supercooled layers colder than — 22°C was also observed, but was much less pronounced than 
at warmer temperatures: we suggest that in cold clouds it is the small droplets in the distribution that freeze into plates and produce 
specular reflection, whilst larger droplets freeze into complex poly crystals. 

The lidar Doppler measurements show that typical fall speeds for the oriented crystals are ^ 0.3ms~^, with a weak temperature 
correlation; the corresponding Reynolds number is Re 10, in agreement with light-pillar measurements. Coincident Doppler 
radar observations show no correlation between the specular enhancement and eddy dissipation rate, indicating that turbulence does 
not control crystal orientation in these clouds. 
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1 Introduction 

It has been recognised for some time that single pris- 
tine ice crystals often fall with their major axes horizon- 
tal over a certain range of Reynolds number. Polarimet- 
ric radar measurements and in-situ sampling led Hogan 
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et al. (2003a) and Field et ah (2004) to suggest that con- 
centrations of these pristine crystals were often associated 
with the presence of mixed-phase cloud layers, in which 
the crystals were nucleated and grew rapidly by vapour 
deposition. These supercooled layer clouds are a fre- 
quent occurrence in the atmosphere (Hogan et al. 2003b, 
2004) but their microphysics is little studied (Fleishauer 
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et al. 2002), and they are poorly simulated in numerical 
weather models at present (Bodas-Salcedo et al. 2008). 

Vertical lidar measurements are an alternative and 
powerful means of observing these pristine crystal pop- 
ulations. If the crystals are planar (and given the lack 
of active ice nuclei at temperatures warmer than — 8°C 
this is the case for many ice-producing supercooled lay- 
ers, Pruppacher and Klett 1997) their basal facets will 
tend to orient in the horizontal plane, and strong 'mirror- 
like' specular reflections are observed as the lidar beam 
is reflected directly back to the receiver (Piatt 1978, Piatt 
et al. 1978). Simultaneous measurements at vertical and 
at a few degrees from zenith allows reliable identifica- 
tion of oriented crystal populations from the large differ- 
ence in backscatter observed at the two angles (Thomas 
et al. 1990). In this paper we apply this technique to iden- 
tify layers of planar crystals, and investigate some of their 
characteristics using a vertically pointing Doppler lidar at 
the Chilbolton Observatory in Hampshire. This instrument 
has operated continuously since September 2006, and the 
extensive data set allows us to provide a much more thor- 
ough statistical investigation of the conditions in which 
these crystals occur than has previously been possible, as 
well as providing many fascinating case studies. These 
are (to the authors' knowledge) the first measurements of 
Doppler velocities from specularly reflecting crystals, and 
this new information provides valuable insight into the fall 
speeds and likely sizes of the crystals involved. 

The article comprises a review of previous observa- 
tions and theory on crystal orientation and specular reflec- 
tion, followed by details of the instruments used, a case 
study of a supercooled layer cloud precipitating oriented 
pristine crystals, and a statistical analysis of 17 months 
of continuous lidar data. A brief discussion concludes the 
paper. 
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2 Review of previous observations and theory 

Although the orientation of ice crystals as they fall is not 
yet fully understood, we can identify some general fea- 
tures. Very small crystals ~ 10/im in diameter are affected 
by rotational Brownian motion (thermal noise) and have 
an approximately random orientation as they fall. Drag 
is produced as vorticity generated at the crystal surface 
diffuses away (Batchelor 1967), and orienting torques are 
weak (Katz 1998): these crystals are unlikely to produce 
significant specular reflection. As the crystals grow big- 
ger and the Reynolds number Re = vD / increases to 
~ 1 (I^ here is the diameter of the crystal, v its fall speed, 
and Uk the kinematic viscosity of the air), the time taken 
to fall a distance equal to the particle diameter t = D/v 
becomes comparable to the time taken for the vortic- 
ity generated at the crystal surface to diffuse that same 
distance t ^ D'^/uk, and the vorticity is 'left behind' as 
it falls. Stable, symmetrical vortices are formed behind 
falling discs (Willmarth et al. 1964, Pitter et al. 1973); 
numerical experiments by Wang and Ji (1997) confirm 
that these standing eddies also form behind hexagonal 
plate and broad-branch ice crystal shapes. If the crystal 
is inclined at an angle to the horizontal plane the asym- 
metrical drag acts to reorient the plate into the horizontal 
position: this reorientation is stabilised by viscosity, so 
oscillations around the horizontal are damped, and the 
particle falls stably and steadily (Willmarth et al. 1964). 
Measurements of natural ice crystals by Kajikawa (1992) 
confirm that horizontal orientation occurs for a variety of 
planar crystal types (all of the PI and P2 types defined by 
Magono and Lee 1966). As the crystals grow even larger, 
the inertial forces are too strong to be damped, and the hor- 
izontal orientation becomes unstable: this initially affects 
the wake downstream producing a slow pitching motion, 
whilst at larger Reynolds number the vortices begin to 
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shed fluid behind the falhng crystal leading to more irreg- (2003) used integrated backscatter measurements to iden- 
ular tumbling motions (Willmarth et al. 1964, Kajikawa tify regions of strong specular reflection, and found it to 
1992). Circular discs fall steadily between 1 < Re < 100; occur exclusively between -10 and — 20°C for 20 days of 
for the natural ice crystals measured by Kajikawa the max- lidar data. Thomas et al. (1990) stepped between and 3° 
imum Reynolds number for stable orientation was in the from zenith to identify regions of cirrus affected by spec- 
range Re 0^ 50-100, and depended on the crystal habit. ular reflection. Collating statistics from cirrus cases over 



Provided that enough crystals are falling in the pre- 
ferred Reynolds number regime, and they have sufficient 
faceted surface, strong specular reflection is possible. We 
note that even branched and dendritic planar crystals typ- 
ically have a faceted plate at their centre (eg. Bentley 
and Humphreys 1964) which may allow specular reflec- 
tion to occur. Piatt (1977) first observed layers of unusu- 
ally strong backscatter in a mixed-phase altostratus cloud. 
Model calculations (Piatt 1978) for populations of ice 
discs with various flutter angles showed that specular 
reflection was able to explain the vertical lidar obser- 
vations. To clinch the interpretation, Piatt et al. (1978) 
made scanning observations in a thin mid-level ice cloud 
where very strong backscatter was observed at vertical 
with almost almost no depolarisation (specular reflection 
at normal incidence does not rotate the plane of polarised 
light); scanning the lidar 2° away from vertical decreased 
the measured backscatter by an order of magnitude, and 
increased the depolarisation ratio to 25%. At the same 
time, a sub-sun (similar to a light pillar, see below) was 
reported by a research aircraft flying above, confirming 
the presence of specularly reflecting crystals. 



3 years they found specular reflection was present approx- 
imately half the time at all heights. They also scanned 
in Imrad intervals around zenith in one case study, and 
inferred from the fall-off in backscatter that the angular 
wobble of the crystals was approximately 0.3°. Young 
et al. (2000) observed 8 mid-level mixed-phase clouds and 
found specular reflection in 20% of layers between —10 
and — 20°C. Sassen and Benson (2001) analysed long term 
statistics of depolarisation ratio in ice cloud, and found 
a minimum in the range —10 to — 20°C which was not 
present in their analysis of off-zenith measurements, indi- 
cating the frequent occurrence of specular reflection at 
these temperatures. Noel and Sassen (2005) used scanning 
lidar measurements of backscatter and depolarisation and 
compared them to forward modelled values to estimate the 
likely range of crystal flutter; their retrieved angular dis- 
tribution half- widths were typically 1-2°. From space Hu 
(2007) analysed returns from the CALIPSO lidar which 
initially pointed 0.3° from nadir, and estimated that half 
of all optically thick ice cloud was affected by specular 
reflection, with important repercussions for retrievals of 
cloud optical properties. The statistics from these studies 
are summarised in table I. 



Since then, a number of similar observations have 
been made. Sassen (1984) observed strong backscatter and 
low depolarisation in the ice falling beneath a supercooled 

liquid layer embedded in a deep orographic cloud. Piatt A number of theoretical calculations have been per- 
et al. (1987) identified layers of ice cloud with anoma- formed to estimate the specular backscatter from a popu- 
lously large backscatter and low depolarisation (integrated lation of oriented crystals (Piatt 1978, Popov and Shefer 
through the layer). Similarly, Hogan and Illingworth 1994, Iwasaki and Okamoto 2001); however knowledge 
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Table I. Previous statistics of specular reflection 



Study Sample Results 

Piatt et al (1987) 48 mid-latitude cirrus cases 42% of profiles with -20 > Tm > -30°C affected, 

20% of profiles with -30 > Tm > -40°C 
+11 tropical cases Trop. cases: no specular reflection observed. 

(nb. Tm = mid-cloud temperature) 
Thomas et al. (1990) Accumulated mid-latitude Present in approximately 50% of cloud at all heights 

cirrus cases from 3 years most of the observations clustered around 9 km altitude 

Young et al. (2000) 8 mid-latitude case studies Observed 20% of time in layers -10 > Tm > -20°C, 

mid-level, mixed-phase did not occur at other temperatures 

Hogan and Illingworth (2003) 20 days continuous data Observed 18% of time in cloud layers at 4 km, 

not observed at cirrus altitudes 
Noel and Sassen (2005) 6 cirrus cases Observed 22% of time for T > -20°C, 

35% in cold clouds T < -30°C 
Hu (2007) 1 month of CALIPSO data Half of all optically thick ice cloud profiles affected 

Chepfer et al. (1999) 3 months POLDER data 37-50% of all ice-only cloud show specular ghnt 

Noel and Chepfer et al. (2004) POLDER, 31 ice cloud cases 80% of ice-only cloud showed specular glint 



of the amount of faceted surface on the crystal, and distri- 
bution of crystal orientations is required, and these quan- 
tities are not well constrained at present. As an example, 
to produce a backscatter of 10~^m~^sr~^ only requires a 
concentration of ~ lm~^ ice discs 500/im in diameter if 
they are all perfectly aligned. However in practice a dis- 
tribution of crystal orientations must always exist, and the 
above observations show that it will be much wider than 
the lidar field of view (66/irad for our monostatic Doppler 
lidar). This means that the true planar crystal population 
is much larger (1000 x larger for a uniform 2° flutter), and 
only a small fraction of the crystals are reflecting the lidar 
beam back to the receiver at any one time. Likewise, if 
only part of the crystal surface is faceted, a correspond- 
ingly larger crystal population is required. Linking crystal 
concentrations to observed backscatter values is therefore 
very difficult. 

In addition to lidar measurements, passive measure- 
ments of radiances from space can also be affected as 
sunlight glints off oriented crystals: Chepfer et al. (1999) 
analysed 30 days of POLDER data and found that this 
happened in at least 40% of ice cloud. Note this instru- 
ment has a much wider field of view than lidar (6.7km 
pixel size). Noel and Chepfer (2004) looked at 31 ice 
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cloud cases, and found 80% were affected by a specular 
reflection. Looking at the change in radiance around the 
glint direction and comparing to simulated reflectances, 
they estimated around 10% of crystals were horizontally 
oriented on average. Most recently, Breon and Dubrulle 
(2004) used a crystal orientation model, combined with 
radiative transfer calculations and estimated that the frac- 
tion of oriented crystals was typically only 1%. Clouds 
containing supercooled liquid were excluded from these 
studies by detection of the rainbow scattering pattern from 
spherical liquid droplets. 

Finally, Sassen (1980) took photographs of light pil- 
lars produced by artificial lights during snowfall. These 
optical phenomena are produced by oriented crystals 
between the light source and observer, which reflect the 
light from their horizontal facets, producing a pillar of 
light apparently extending vertically from the light source. 
Any crystal wobble produces a spreading out of the pillar, 
as rays of light out of the plane formed by the light pillar, 
the light source and the observer are also reflected (to the 
observer), leading to a finite pillar width. By measuring 
this pillar width, Sassen estimated the characteristic wob- 
ble of the crystals in the snowfall, with observed values 
in the range 0.5-3°. Correlating this width with the mean 

QJRMS 00: 1-26 (2009) 



SPECULAR REFLECTION FROM ICE CRYSTALS 



5 



Table 11. Specifications of Doppler Lidar parameters with a 30 second integration time, and has a 



Wavelength 1.5 /im gate length of 30 metres. This instrument was operated 

Range resolution 36 m 

Integration time 30 s 45° elevation for the case study in section 4, and at 

Maximum range 10 km . , ^ i ... • ^ rr^i 

Maximum Doppler velocity 10 m/s ^^^^ti^^l Statistics m section 5. The 94 GHz radar 

Doppler velocity noise 0.05 m/s (SNR dependent) i i • i i-^ • 

Pulse duration 200 ns operates on a case study basis to preserve the liretime or 

PRE 15 kHz 

Beamwidth 60 mm (at source) 

Divergence 33 /xrad 



Antenna monostatic optic-fibre coupled 



the tube, measuring reflectivity, and Doppler parameters 
at 1 second resolution, with a gate length of 60 metres. 
Full Doppler spectra are also recorded once per minute by 
size of crystals (mostly plate and branched planar crys- ^oth radars when cloud is present; the integration time for 
tals) collected at the ground, and using empirical velocity- these spectra is 1 second, 
diameter relationships, his data showed a minimum in the 



pillar width (and therefore crystal wobble) at a Reynolds 
number of Re 8. 

3 Instruments and method 



A key parameter in this study is the lidar 'colour 
ratio' which we define as the ratio of the lidar backscatter 
measured by the 905nm ceilometer to that measured by 
the 1.5/im Doppler lidar: 



'905nm 



^1.5/2m 



(1) 



The key instrument in this study is a newly developed Colour ratio [dB] — lOlog^o 
1.5/Lim Doppler lidar (HALO Photonics Ltd., Malvern, 

UK). The instrument measures profiles of backscatter and Under conditions where there is no specular reflection we 

mean Doppler velocity resolved over 36 metre range gates expect this ratio to be greater than OdB for ice particles, 

and recorded every 30 seconds to a range of 10 km. ^in^e the imaginary component of the refractive index of 

The lidar has operated continuously since September i^e is 3 orders of magnitude larger at 1.5/xm compared 

2006, and points directly at vertical. Specifications are ^ ^O^nm (n" = 4.7 x IQ-^ vs. n" = 4.3 x 10"^ Warren 

given in table II. In addition to this instrument, a 905nm 1984). This results in significant absorption of photons as 

lidar ceilometer records backscatter profiles at approxi- they travel through the ice crystal, and a corresponding 

mately the same range and temporal resolution (Illing- reduction in the backscattert. As a guide to the magnitude 

worth et al. 2007), and points 4° off vertical. This lidar is this effect, consider a ray of light as it undergoes a series 

calibrated using measurements of the integrated backseat- internal reflections through the ice crystal and back to 

ter in optically thick stratocumulus cloud as described by lidar telescope. According to the Beer-Lambert rule, 

O'Connor et al. (2004). Mie calculations for liquid water intensity will be reduced by a factor exp(-47rn"x/A), 
droplets indicate that the 1.5/im lidar can be calibrated in 

a similar way: details of this are given in Appendix A. ^We note that the attenuation of the hdar beam as it passes through a 

unit volume of cloud depends only the extinction of the ice particles 

Radar measurements were made usins 35- and 94- within the cloud, which should be the same at both wavelengths (equal to 

twice the projected area of the crystals per units volume of air, assuming 

^TTT^ 1 A 4- A -uju scomctric optics). As a result the ratio of the measured attenuated 

GHz Doppler radars: both mstruments are described by ^^ekscatters should be identical to the ratio of the true (unattenuated) 

1 -1 backscatters, provided that the ice cloud sampled by the two lidar beams 

Illmgworth et al (2007). The 35GHz radar provides is the same (ie. it is horizontally homogeneous on scales of a few 

hundred metres). We have neglected attenuation by aerosol particles, 

reflectivity (at two orthogonal polarisations) and Doppler and multiple scattering between ice crystals. 
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off-zenith data 
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Colour ratio [dB] 

Figure 1 . Distribution of colour ratios (905nm / 1 . 5/im backscatter) 
collected in ice cloud while both lidars were pointing off-zenith. 
Frequencies are normalised such that the total area under the curve 

=10. 

SO for a total path length ofx = 200/im, the backscat- 
ter at A = 1.5/im is reduced by 3.4dB, whilst light at the 
ceilometer wavelength is essentially unaffected. We con- 
clude that this differential absorption is likely to be signifi- 
cant for ice crystals in the atmosphere which are typically 
a few hundred microns to a few millimetres in size, and 
expect positive colour ratios > OdB. This effect could in 
principle be used to measure crystal size in ice clouds, and 
will be explored in a separate publication. 

To check the analysis above is correct, we have col- 
lected 9 days of backscatter data with both lidars pointing 
off-zenith. During this period (25 February to 4 March 
2008) various types of ice cloud were observed includ- 
ing cirrus, altocumulus with ice fall streaks below, and 
thick frontal ice cloud. Quicklooks of the lidar data may be 
viewed at www.met.reading.ac.uk/radar/realtime/. From 
this data figure 1 shows the histogram of the measured 
colour ratio in the ice-phase regions. It is apparent that 
the absorption effect is significant, with 95% of the data 
corresponding to colour ratios between +3 and +12 dB 
(ie. backscatter at 1.5/im reduced by a factor 2 to 16). 
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Almost no data were observed with colour ratios < OdB 
— the 1% of pixels that were, revealed themselves on 
closer inspection to be isolated points resulting from occa- 
sional artefacts recorded by the ceilometer. 

These off-vertical measurements contrast sharply 
with our observations with the 1.5/im Doppler lidar point- 
ing directly vertical. In such cases we frequently observe 
colour ratios < dB, in some cases as low as — 25 dB (ie. 
300 X more backscatter at vertical). We conclude that these 
negative colour ratios must correspond to cases where the 
ice crystals are highly aligned in the horizontal, ie. the 
specular reflection from the horizontal crystal facets pro- 
duces a strong signal in the vertical Doppler lidar data 
which is not present in the off-vertical ceilometer data, 
leading to the negative colour ratios. In what follows we 
will use this condition to identify regions of cloud where 
specular reflection is occurring (negative colour ratios) 
and where it is not (positive colour ratios). 

4 Case study, 18 May 2008 

4. 1 Vertical measurements 

Here we present an example of a persistent mid-level 
mixed-phase cloud associated with a weak occlusion mov- 
ing southwest across England on the 17 and 18 May 
2008. Patchy mid-level cloud cover was first observed 
over Chilbolton by the vertically pointing 94GHz cloud 
radar at approximately 06 UTC on the 17*^; by 15 UTC 
this had become a well defined layer, with cloud top at 
4000m altitude. At the same time low-level stratocumu- 
lus was obscuring the mid-level cloud from view of the 
lidars; in addition some of the ice fall streaks appeared 
to be seeding the stratocumulus, producing drizzle. How- 
ever, by early morning on the 18 May the low cloud began 
to break up allowing an unobscured view of the mixed- 
phase cloud above: figure 2 shows the time series from the 
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ceilometer and Doppler lidar for the period 01-10 UTC. 
The off-zenith ceilometer clearly picks out the location 
of the thin supercooled liquid water layers, with a well 
defined primary liquid layer at the top (3 800-4 100m alti- 
tude), and intermittent secondary liquid layers embedded 
in the l-2km deep ice fall streaks. The concentration of 
ice crystals is much lower than the concentration of cloud 
droplets, and so the backscatter is an order of magnitude 
smaller in the ice virga than in the supercooled liquid 
layer; the fall streak structure of the sedimenting ice crys- 
tals is clearly visible here. In the lowest 1km weak returns 
from the boundary-layer aerosol are visible in the lidar 
time series; the radar meanwhile is insensitive to these tiny 
aerosol particles. 

The OZ sounding from Herstmonceux (125km east of 
Chilbolton) is shown in figure 3, and the water-saturated 
layers at 3 800-4 100m and 2500-2700m seem to be con- 
sistent with the ceilometer imagery. The temperature at 
the top of the uppermost liquid layer was — 15°C; a slight 
temperature inversion was present at cloud top, and the air 
above was much drier. Wind speeds were relatively light; 
directional shear was observed in the ice virga, and is also 
apparent from the shape of the fall streaks in the radar time 
series, figure 2). 

Data from the vertically-pointing Doppler lidar are 
also shown in figure 2, and it is immediately apparent 
that the backscatter in the ice virga is typically one order 
of magnitude (10 dB) larger than that recorded by the 
ceilometer, in spite of the increased absorption at 1.5/im. 
This convinces us that oriented planar crystals are present 
in significant numbers in the virga. These crystals must 
be nucleated in the primary supercooled layer, growing by 
deposition in the vapour-rich conditions. The temperature 
in this water- saturated layer was between —12.5 and 
— 15°C leading to a large supersaturation over ice (see 
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Figure 3. Herstmonceux ascent at UTC on 18 May 2008. Panels 
clockwise from top left show: air temperature (black) and dew- 
point (grey); relative humidity with respect to liquid (grey) and ice 
(black); wind speed; wind direction in degrees. 



figure 3); Ryan et al (1976) found skeletal plates (at 
— 13) and stellar crystals (at —15) grew rapidly in such 
conditions. The aspect ratio of their crystals was measured 
to be more than 10: 1 by the time they had reached 100/im 
in size, and given the discussion above this should lead to 
a strong reorienting torque (indeed Ryan et al noted that 
nearly all of the crystals settled flat-on as they fell onto the 
sampling slide). Because of the large vapour excess, the 
crystals do not grow in such a controlled way as simple 
hexagonal plates; however there is still likely to be ample 
faceted surface for specular reflection to occur from. 

Unfortunately the microwave radiometers at 
Chilbolton were not operating on this day. However 
Korolev et al. (2007) have developed a correlation 
between the depth of mid-latitude supercooled layers 
and their liquid water path from aircraft data, and from 
this we estimate a value of only 40gm~^; this low 
value is consistent with previous measurements of similar 
supercooled clouds (eg. Hogan et al. 2003 a,b). As a 
result riming is likely to be very light or not at all, and 
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Ceilometer backscatter, off zenith 
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Figure 2. Lidar and radar observations, morning of 18 May 2008. Top panel shows the backscatter measured by the ceilometer pointing at 
4° from zenith; second panel shows the backscatter measured by the Doppler lidar pointing directly at zenith. Colour scale is in logarithmic 
units ie. logiQ(m~^sr~^). Third panel shows the vertical velocity measured by the Doppler lidar (negative velocities correspond to 
downward motion). Final two panels shows the 94-GHz radar reflectivity and vertical Doppler velocity. 
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the observed strong specular reflection supports this 
assertion. 

In contrast to the hdars, the vertically pointing 
94GHz radar (figure 2) is insensitive to the supercooled 
liquid droplets in the mixed-phase layer, and is domi- 
nated by returns from the much larger ice crystals; again 
fall streaks are clearly visible. The reflectivity of the ice 
crystals falling from the primary supercooled layer is 
observed to vary in the range -15 to dBZ, suggesting ice 
water contents of ^ 0.01-0. lgm~^ (Hogan et al 2006); 
this range is consistent with in- situ observations of mid- 
latitude mixed-phase layer clouds (Carey et al 2008). 

Figure 2 shows the vertical velocity from the liquid 
droplets and ice particles as measured by the Doppler 
lidar. We have calculated the average Doppler velocity in 
the uppermost liquid layer for the period shown in figure 
2: this is approximately +0.02ms~^, indicating that any 
large scale ascent is very gentle. However it is clear that 
there are smaller scale vertical air motions present - at 
cloud top there is convective overturning with peak-to- 
peak vertical motions of 0.5ms ~^ as the cloud cools 
radiatively to space. Photographs of the sky in figure 4 
confirm the cellular structure. Waves (approximately 10 
minute period) are also present - these are most visible in 
the earliest part of the time series, and also in the vertical 
radar Doppler velocity, also shown in figure 2. 

To get an estimate of the terminal velocities of the 
oriented ice crystals, we have isolated areas with colour 
ratio < OdB and calculated the mean lidar Doppler veloc- 
ity over the 9 hour period. By averaging over several hours 
we expect to remove the mesoscale air motions observed 
in the detailed time series and obtain approximate values 
for the still- air velocities of the crystals. We estimate this 
to be -0.23ms"^ Sassen's (1980) light pillar data indi- 
cated that there is likely to be a Reynolds number (Re ^ 8) 
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Table III. Estimated Reynolds number and crystal diameter corre- 
sponding to fall velocity of 0.23ms~^ 

Assumed habit Re D [/xm] 

Hexagonal Plates 4 400 

Sector Plates 5 550 

Broad-branch 5 550 

Plates w/extensions 6 650 

Dendritic 9 1000 

Stellar 11 1200 



where the crystal flutter is minimised, and we expect 
that these crystals will provide the strongest weighting 
to the Doppler velocity measurements. Heymsfield and 
Kajikawa (1987) measured the fall velocity of planar ice 
crystals collected at the ground in northern Japan (1 km 
above sea level), and formulated velocity-diameter rela- 
tionships from those observations. Using the pressure cor- 
rection method of Beard (1980) we have used this data to 
calculate the ranges of Reynolds number and crystal diam- 
eter corresponding to a fall speed of 0.23ms~^ for various 
assumed crystal habits: these values are tabulated in (III). 
We find Re ^ 4, D ^ 400/im for the most compact habit 
(hexagonal plates) increasing to Re 11, 1) 1200/im 
for the most open (stellar crystals). This range of Reynolds 
numbers is broadly consistent with Sassen's observations, 
but may be more robust since we are measuring velocities 
directly rather than median crystal diameters at the ground 
(this may be significant since the distribution of crystal 
orientations is not the same for all crystal diameters). Our 
observations are also consistent with Kajikawa's (1992) 
observations since the Reynolds numbers estimated here 
are smaller than the critical values where the fall motion 
becomes unstable. 

Figure 5 shows the linearly averaged lidar backscat- 
ter profile from the ceilometer and Doppler lidar between 
0530 and 0545 UTC (a period when the cloud appears 
to be relatively steady- state and horizontally homoge- 
neous). The profiles reveal that that specular reflection is 
maximised at approximately 3000-3300m (colour ratio 
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Figure 4. Photographs of the sky every 30 minutes from 0420 to 0650 UTC 
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Figure 5. Averaged hdar backscatter profiles for the period 0530- 
0545 UTC. Grey line is the profile from the off-zenith ceilometer, 
black line is the profile from the vertically pointing Doppler lidar. 
Shading indicates the difference in backscatter measured by the two 
instruments: dark grey shading corresponds to specular reflection 
from oriented crystals; light grey shading corresponds to evaporat- 
ing crystals which have lost their faceted surfaces. 



= -15dB), decreasing steadily at lower altitudes. The 
backscatter profiles in the supercooled liquid are well 
matched, with slightly higher returns from the ceilome- 
ter due to the difference in backscatter to extinction ratio 
between the two wavelengths (see appendix A). This is 
evidence that the lidars are well calibrated, and that the 
attenuation through the ice virga is the same for both lidar 
beams, despite the enormous difference in backscatter. 

The regions near cloud base with positive colour 
ratios +4 dB) are noted with interest. The fall off in 
lidar and radar backscatter in this region indicates that it is 
subsaturated with respect to ice. In these conditions evap- 
oration tends to proceed from crystal edges since this is 



where the molecules are most weakly bound. If the sub- 
saturation is very slight, evaporation occurs slowly in dis- 
crete layers, allowing the crystal to maintain its faceted 
structure (Hallett et al. 2001). If the crystal falls into 
drier air however, multiple steps form at the edges (Mason 
1971) and the evaporation rate becomes limited by the dif- 
fusion of the molecules through the air, leading to strong 
vapour gradients at the edges and comers (Nelson 1998, 
Westbrook et al. 2008), producing smooth, rounded crys- 
tals (Nelson 1998). Rounding also appears to be more 
favourable at warmer temperatures (Hallett et al. 2001). 
Given these results it seems clear that the regions of pos- 
itive colour ratios are cloud where evaporating crystals 
have lost their faceted surfaces and become rounded, so 
specular reflection can no longer occur^ . This is consis- 
tent with calculations by Mishchenko et al. (1997) who 
showed that flat-faced discs could produce strong spec- 
ular reflection, whilst thin oblate spheroids (ie. slightly 
rounded faces) could not. 

An alternative explanation might be that the crystals 
have evaporated to the point where the crystals are too 
small to be horizontally oriented (tens of microns in 
diameter, see section 2). However, given that the radar and 
lidar Doppler velocities in this region are consistent with 
crystals several hundred microns in size, we discount this 
idea. 



^This idea also appears to explain the vertical lidar observations of two 
mid-level clouds by Uchino et al. (1988) who were puzzled when they 
found that that whilst the returns from the middle of their clouds were 
characteristic of specular reflection (very high backscatter, very low 
depolarisation ratios), the returns from the base of the virga were char- 
acteristic of normal ice cloud (much lower backscatter, depolarisation 
ratios increased to ^ 30%), and that this was co-located with a rapid 
decrease in relative humidity with respect to ice. 
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Identification of these regions where specular reflec- 4ioo 
tion is 'switched off may therefore be a useful observa- 
tional tool to determine if the crystals are evaporating in 
layers or losing their shape. The latter case is also evi- 
dence that the whole crystal surface is close to ice satura- 
tion (Nelson 1998), and that the crystal evaporation rates S: 
are governed by the standard electrostatic approximation 
(Mason 1971, Westbrook et al 2008). 



To investigate the initial growth of the ice at the top 
the cloud we have averaged the vertical 94-GHz radar 
Doppler velocity and reflectivity over the whole period 
(01-10 UTC) for the top 600m, as shown in figure 6. At 
this range of altitudes the radar data suggest that there 
is a relatively consistent production of ice. The averaged 
profiles show that the reflectivity and Doppler velocity of 
the ice particles increase steadily in magitude through the 
supercooled layer (the top 300m) reaching values of ^ 
— 3dBZ and — 0.45ms~^ respectively. Note that the radar 
velocities are larger than measured by the Doppler lidar 
in the specular ice, reflecting the weighting of the radar 
towards the heaviest ice particles. The crystal diameter 
corresponding to this velocity is rather sensitive to the 
assumed habit, ranging from 800/im for hexagonal plates 
to 3500/im assuming stellars (Heymsfield and Kajikawa 
1987). We note that over the depth of the supercooled 
layer the reflectivity increases by about 6dB, which corre- 
sponds to a doubling of the average crystal mass over that 
distance (since reflectivity is, to first order, proportional to 
mass^) whilst Ryan et al. (1976) observed crystal masses 
to increase by a factor of ten between 50 and 150 seconds 
of growth. The comparison is not so clear-cut because of 
non-Rayleigh scattering effects and sedimentation, but the 
implication is that the crystals are not simply nucleated at 
cloud top and grow whilst steadily falling falling through 
the water saturated layer, but are also mixed through the 
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Figure 6. Averaged radar reflectivity and velocity profiles for the 
period 01-10 UTC for the top 900m of cloud. 



layer as the air is gently turned over, smoothing out some 
of the reflectivity structure. It may be that crystals are also 
nucleated throughout the depth of the liquid layer rather 
than only at the top, which would have a similar effect. 
Below 3800m there appears to be relatively little change 
in reflectivity or velocity, suggesting little growth is occur- 
ing - presumably the relative humidity falls off sharply to 
near ice saturation here (there is some evidence for this in 
the radiosonde profile, figure 3). 

Unlike Shupe et al (2004) we were unable to detect a 
signature from the supercooled liquid droplets in the radar 
Doppler spectra (not shown): a single mode was observed, 
corresponding to the sedimenting ice crystals. Most likely 
this is because the liquid water content in this cloud is too 
low to have a detectable radar signature: the case studied 
by Shupe et al had a peak liquid water path of 500gm~^ 
in a 1km thick layer with droplets ascending at lms~^; 
our example is much thinner and more weakly forced, and 
our estimated liquid water path is an order of magnitude 
smaller. 

Figure 7 shows the turbulent eddy dissipation rate e, 
estimated from the vertically pointing radar as described 
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by Bouniol et al. (2003). The standard deviation of 1- 
s snapshots of mean Doppler velocity over 30s periods 
are combined with wind speed estimates (from the Met 
Office 12km model forecast), and e is inferred from the 
variability in the vertical velocity over the 30s sample 
time. In most of the virga the air appears to be quite 
still with dissipation rates of order 10~^ - 10~^m^s~^, 
including the evaporation zone at the base of the virga; 
however in the top 500m of the cloud, this peaks at 
e^3xl0~^m^s~^, presumably reflecting the radiative 
cooling of the liquid droplets and associated overturning. 
This turbulent mixing seems to extend down into the 
top few hundred metres of the virga; despite this the 
lidar data indicates that strong specular reflection is still 
occurring in this region. This indicates that turbulence is 
unlikely to be a strong influence on crystal orientation, 
and this is consistent with previous theoretical analysis 
(Klett 1995): for e = 3 x IQ-'^m^s-^ the Kolmogorov 
microscale (the size of the smallest eddies) is 1.8mm, 
and the corresponding velocity perturbation across a 1mm 
ice crystal is of order 0.005ms~^, much smaller than its 
fall speed. Interestingly, as well as the turbulence in the 
liquid layer at the top, there are also regions of enhanced 
e co-located with the more intermittent supercooled layers 
embedded in the ice virga (eg. 0630-0700 UTC, 3200m), 
perhaps associated with the flow of latent heat as these 
layers are seeded by the crystals falling from above. 

4.2 Polarimetric radar measurements at 45° 

In addition to the vertical and near- vertical radar and lidar 
data shown in figure 2, we also made simultaneous polari- 
metric radar measurements with the 3 5 GHz cloud radar 
dwelling at 45° (due West; the cloud is advecting in from 
the North East). The corresponding time series is shown 
in figure 8. The radar reflectivity structure is very simi- 
lar to that recorded by the vertical 94GHz radar, giving 
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Figure 7. Eddy dissipation rate e estimated from the standard 
deviation of 1-s mean radar Doppler velocities over 30s periods. 
Grey scale is in logarithmic units, ie. log^Q(m^s~^). 

US confidence that the cloud is essentially homogeneous 
within a 4km range of Chilbolton, and that we are sam- 
pling similar ice particles. Also shown in this figure is the 
differential reflectivity Zdr. This parameter is the ratio 
of the co-polar radar reflectivity for horizontally polarised 
waves vs. those polarised at 45° from horizontal (and 
perpendicular to the direction of propagation). Calcula- 
tions (Appendix B) show that thin pristine crystals which 
are horizontally oriented produce a strong Z^r signature 
(up to 3.5 dB), whilst we expect randomly oriented crys- 
tals and low-density/irregular particles have Z^r ?^ dB. 
Crucially, if there are only a few pristine crystals mixed in 
with large irregular particles, the Z^r signature is masked 
(see Hogan et al. 2002,2003a; Eield et al. 2004), because 
the reflectivity at both polarisations is dominated by the 
large irregular particles. 

The observed values of Zdr demonstrate that the 
15 dB of specular reflection seen in the lidar data is not 
the result of a few perfectly formed and oriented crys- 
tals, but reflects a broader population of planar crystals 
falling close to horizontal. Most of these are probably 
imperfect as a result of growing at water- saturation, and 
have a distribution of orientations centred around hori- 
zontal; but collectively these crystals have enough faceted 
surface normal to the lidar beam to produce a strong 
specular signal. Since crystals grown in this temperature 
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Figure 8. 45° radar observations. Top panel shows the reflectivity recorded by the 35-GHz cloud radar; lower panel shows the differential 
reflectivity Zdr- High values of Zdr correspond to returns from strongly oriented planar crystals. 



3.5 



2.5 



DC 
Q 

N 1.5 



0.5 





1 1 


1 1 1 -pr 















... <o 










1 1 





0.2 



0.4 0.6 
Volume fraction 



0.8 



Figure 9. Zdr for horizontally oriented planar crystals observed at 
45° elevation as function of volume fraction. Lines indicate range 
of crystal aspect ratios observed by Ryan et ah 1976 (0.02-0.14 
in steps of 0.04 from darkest to lightest). Typical ice/air volume 
fractions for some crystal types are indicated, inferred from Ryan 
et al and Heymsfield and Kajikawa (1987). 



range are typically very thin, Zdr is primarily sensitive to 
the crystal ice/air volume fraction (the volume of the ice 



crystal divided by a spheroid of the same overall dimen- 
sions), and calculated values for perfectly horizontal crys- 
tals are shown in figure 9 . The detail of this calculation 
is described in Appendix B; we have marked on the parts 
of the curve corresponding to hexagonal plates and small 
stellar crystals as inferred from Ryan et a/.'s (1976) data. 
Given the cloud top temperature we anticipate growth ini- 
tially of simple plates, evolving into branched stellar crys- 
tals after 1-2 minutes of growth, implying Zdr ^ 2-3 dB; 
the observed differential reflectivity measurements seem 
to back up this expectation, with much of the virga hav- 
ing values in this range. There is a fascinating amount 
of structure in the Zdr observations however; presum- 
ably this is largely the result of size sorting in the fall 
streaks yielding streaks of high Zdr where the crystals 
are small and dense, and lower Zdr where the crystals 
are larger and more tenuous. The fall streak structure is 
similar to that observed in supercooled stratocumulus by 
Field et al. (2004). We also note that the exact crystal habit 
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is quite sensitive to the growth temperature, and given the 
range of temperatures spanned by the supercooled layer, 
it is quite likely a range of planar habits develop simulta- 
neously, and are then sorted as they fall. The very largest 
values of Z^r 3dB) seem to be found only close to the 
very top of the cloud: presumably these are small plate- 
like crystals in the early stages of growth, with a volume 
fraction close to unity (see figure 9); the crystals lower 
down will likely have been growing for longer, developing 
tenuous branches, perhaps aggregating, and will therefore 
will have a lower volume fraction and lower Zdr- The 
observation of strong differential reflectivity even at the 
very base of the virga where the crystals are evaporat- 
ing is further evidence that the lack of specular reflection 
in the evaporation zone is due to the crystals becoming 
rounded (whilst still maintaining their horizontal orien- 
tation) rather than their orientation being disrupted, for 
example by turbulence associated with evaporating cool- 
ing. 

5 17 months of statistics 

We now present a statistical analysis of specular reflection 
from 17 months of continuous observations. A number 
of the studies cited in section 2 have calculated statistics 
on the occurrence of specular reflection, as summarised 
in table I. Many of the results are inconsistent with one 
another, and all but Young et al. (2000) and Hogan and 
Illing worth (2003) have effectively excluded profiles con- 
taining supercooled liquid. Unfortunately these latter two 
studies sampled quite a small number of clouds. It is there- 
fore desirable to calculate some statistics of occurrence 
for specular reflection using our long continuous data set, 
and since we are able to range-resolve the position of the 
oriented crystals we can readily investigate their occur- 
rence below supercooled layers, as well as in ice- only 
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profiles. We note that our results are distinct from those 
derived from POLDER since our 30s x 36m sample vol- 
ume is much smaller. The colour ratio and Doppler mea- 
surements allow us to explore more deeply some of the 
properties of these oriented ice crystals, such as how fast 
they fall, and what part of the distribution they likely rep- 
resent. 

The data were analysed as follows. Boundary layer 
aerosol and precipitation were excluded by removing 
targets which produced a continuous profile of backscatter 
down to the lowest range gate. This effectively excludes 
supercooled stratocumulus from our analysis, and may 
also exclude occasional mid-level cloud where ice fall 
streaks extend into the boundary layer aerosol. In general 
however this procedure is quite effective at isolating non- 
precipitating cloud. As an additional condition we restrict 
our attention to cloud where the temperature is colder 
than — 2.5°C in order to exclude melting ice particles; 
the temperature is estimated from the Met Office's 12km 
forecast model output over Chilbolton (see Illingworth 
etal 2007). 

We identify cloud as range gates where the backscat- 
ter is greater than 5 x 10~^m~^sr~^ from both lidars. 
The presence of supercooled liquid droplets is inferred 
by backscatter from the off- vertical ceilometer exceeding 
10~^m~^sr~^ and this range gate is flagged and removed. 
The backscatter from the three range gates below this are 
also removed from the analysis since these may also con- 
tain liquid, and all the data from range gates above is 
also discarded since liquid cloud may well have signifi- 
cant structure on the scale of hundreds of metres and as a 
result there may be some differential attenuation between 
the zenith and off-zenith lidar beams, which could per- 
haps lead to an ambiguous interpretation of colour ratios 
measured from ice cloud above. The remaining ice cloud 
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Figure 10. Percentage of ice cloud which is identified as containing 
specularly reflecting ice crystals, binned as a function of tempera- 
ture. A distinct peak is apparent close to — 15°C, where 40% of ice 
cloud measured by both lidars is observed to be influenced by spec- 
ular reflection. At temperatures colder than — 30°C this fraction is 
only 10%. These figures are hkely underestimates, and the dashed 
line shows the fraction of cloud affected if all colour ratios < +3dB 
are counted (see text). The total number of ice pixels in each tem- 
perature bin ranged from 3 x 10^ in the coldest bin, to 1.5 x 10^ 
in the warmest. 



is divided into data with colour ratios < dB (specular 
reflection) and data with colour ratios > dB ( 'normal' 
ice cloud). Note that the lidar signal is attenuated even in 
glaciated clouds, and whilst it is often possible to obtain a 
complete backscatter profile through a relatively shallow 
cirrus cloud or altocumulus virga, for deep stratiform ice 
clouds we are typically only able to gather statistics from 
the lowest kilometre or two of cloud. Convective clouds 
are effectively excluded since the lidar is not able to see 
into the liquid cloud base. This analysis was performed for 
a total of 503 days (approximately 17 months) of contin- 
uous lidar observations between 12 September 2006 and 
31 January 2008. 



5 . 1 Occurrence by temperature 

Using the method described above we have calculated 
the fraction of cloud that is observed by both lidars and 
which is identified as containing specularly reflecting ice 



crystals. Overall we find that of the 1.55 x 10^ ice cloud 
pixels analysed, 26% were affected by specular reflection. 
Breaking this fraction down as a function of temperature 
reveals further insight, as shown by the horizontal bars 
in figure 10 - it is apparent that most of the specular 
reflection is concentrated at temperatures warmer than 
— 25°C, with a maximum in the region of — 15°C where 
over 40% of ice cloud observed is affected by specular 
reflection. These observations seem consistent with the 
expectation that planar monocrystals typically grow in 
the temperature range —9 to — 22°C (Pruppacher and 
Klett 1997). There is no obvious 'gap' corresponding to 
the window of column/needle growth between —3 and 
— 8°C, but since there are few ice nuclei at these warm 
temperatures anyway, the specular in this region is likely 
planar crystals formed at colder temperatures sedimenting 
in from above. The maximum occurrence around — 15°C 
seems sensible, since thin planar crystals are well known 
to grow in this temperature region. 

At temperatures between —30 and -40°C only 
around 10% of cloud is affected by specular reflection. 
This is a significantly smaller fraction than observed in 
cirrus clouds by Thomas et al. (1990), although their scan- 
ning lidar arrangement may be a more sensitive tool to 
detect weak specular reflection, since such a setup is not 
complicated by differential absorption effects (see below). 

Figure 11 shows the observed distribution of lidar 
colour ratios as a function of temperature. The lack of 
bimodality around the OdB mark suggests that in fact 
specular reflection is likely to be occurring weakly in 
some of the cloud where we measure positive colour 
ratios. This is a result of the increased absorption at 
1.5/im compared to 905nm as discussed in section 2. 
Specifically, this would correspond to cases where the 
specular reflection component of the backscatter is of 
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comparable magnitude to the backscatter from the other 
ice particles in the cloud, and the overall colour ratio is 
positive (our off-zenith measurements indicated that ice 
particles which are not affected by specular reflection 
yield colour ratios in the range +3 to +12dB). 

The distinction between specular and non- specular 
ice is most obvious at the mid-level temperatures between 
-10 and -20°C, where a long straight tail extending 
out to — 20dB is apparently distinct from a more bell 
shaped distribution for normal ice cloud, with the two 
curves intercepting at approximately +3dB. At tempera- 
tures colder than — 30°C the distribution is much narrower 
with relatively few colour ratios < OdB, and the magni- 
tude of the specular reflection (as characterised by the 
colour ratio) in these cold clouds is much weaker than at 
-15°C. 

We note that if it were assumed that all values less 
than +3dB in fact corresponded to some specular reflec- 
tion occurring in the lidar sample, the overall fraction 
of ice cloud affected would rise from 26% to 40%. The 
occurence by temperature following this assumption is 
shown by the dashed line in figure 10. Our estimates of the 
frequency of specular reflection should therefore be con- 
sidered as lower bounds on the true occurence, and this 
emphasises the large fraction of cloud in which specular 
reflection is possible, especially at the mid-levels. In the 
subsections that follow we will use the dB threshold to 
identify cloud as specular or not, ie. where the effect is 
strong enough to be measureable without ambiguity. 

5.2 Fall speeds 

Applying the method of identification above, we can use 
the Doppler velocity measurements from the vertical lidar 
to study how fast the oriented crystals fall. The frequency 
distribution of Doppler velocities for all cloud affected 
by specular reflection is shown in figure 12, which also 
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Figure 1 1 . PDFs of the lidar colour ratio seperated by temperature. 
At temperatures characteristic of planar growth (especially between 
— 10 and — 20°C) the distribution has a long tail with ratios < OdB 
(specular reflection), whilst at colder temperatures the distribution 
is much narrower and concentrated at colour ratios > OdB (normal 
ice cloud). 

shows the velocity distributions for normal ice cloud and 
liquid droplets. The frequencies are normalised so that the 
area under each curve is equal to 1 . Note that the veloc- 
ity distributions are broadened somewhat by the range of 
vertical air motions in the clouds, as well as by variability 
in the particle terminal velocities. The distribution for liq- 
uid droplets has a mean velocity of 0.00 ms~^, the width 
of the distribution is likely dominated by turbulence of 
the kind observed in figure 7. For normal ice cloud we 
measure a relatively broad distribution of velocities, with 
a mean value of — 0.49ms~^ and standard deviation of 
almost exactly the same magnitude. This suggests that in 
the majority of ice clouds sampled by the lidar, the lidar 
signal is controlled by ice particles of several hundred 
microns in size. In cloud affected by specular reflection, 
the range of crystal fall speeds is rather narrower, with 
a mean value of — 0.32ms~^. This is in keeping with the 
idea that specular reflection is produced by planar crys- 
tals from a relatively limited range of Reynolds number 
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Figure 12. Distribution of Doppler velocity for specularly reflect- 
ing ice crystals (black line), 'normal' ice cloud (solid grey line), 
and supercooled liquid droplets (dashed grey line). Frequencies nor- 
malised so area under the curve = 1. Mean values are —0.32 ms~^, 
— 0.49 ms~^, and 0.00 ms~^ respectively. 



where the fall orientation is stable and the fall speeds are 
relatively slow. The corresponding Reynolds numbers for 
various assumed crystal habits have been calculated using 
Heymsfield and Kajikawa's (1987) relationships. Segre- 
gation of the data by colour ratio (not shown for brevity) 
shows a slight decrease at the most extreme colour ratios 
(10-20%); however this decrease is much smaller than 
the variation between the different assumed crystal habits 
(Re ^ 5 for hexagonal plates vs Re ^ 20 assuming stel- 
lars). In any case it seems that the specular backscatter 
and Doppler velocity is weighted to crystals falling with 
Re ~ 10, with a factor of 2 uncertainty. This is consistent 
with Sassen's light pillar observations, and the case study 
in section 4. 

Plotting the Doppler velocity distribution for specu- 
lar crystals as a function of temperature yields figure 13. 
The average particle fall speed for normal ice cloud was 
calculated by Westbrook and Illingworth (2009) and was 
found to increase systematically with temperature from 
0.3ms"^ at -35°C, to 0.65ms"^ at -5°C, demonstrating 
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Figure 13. Distribution of specular fall speeds as a function of 
temperature, filled contours indicate frequency normalised for each 
temperature bin in 0.2 [sm~^] intervals. White fine with triangles 
indicates the mean values for specular cloud as a function of 
temperature; stars indicate same but for normal ice clouds. 



the influence of larger particles and broader size spec- 
tra at warm temperatures. In contrast to this, the veloc- 
ity distribution for specularly reflecting crystals shows a 
rather weaker variation with temperature, with average 
crystal fall speeds increasing from 0.25ms~^ at — 35°C, 
to 0.3ms~^ at — 15°C; there appears to be a slight broad- 
ening of the distribution for the warmest clouds, with an 
average velocity of 0.4ms~^ at — 5°C. Possibly this corre- 
sponds to planar crystals formed at the mid-levels which 
thicken as they fall into air temperatures characteristic 
of columnar growth, increasing their weight relative to 
their drag, and allowing them to sediment more quickly 
for a given Reynolds number. As in section 4 the mean 
Doppler velocities and implied crystal Reynolds num- 
bers are lower than the maximum stable values found by 
Kajikawa (1992) for all assumed crystal types. 
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5 . 3 Crystal alignment 

The correlation between turbulent air motion and colour 
ratio measurements has also been investigated, to deter- 
mine whether the horizontal orientation of the ice crystals 
is affected by the turbulent eddies. The eddy dissipation 
rate e was estimated as described in section 4, and the 
results were binned by colour ratio and dissipation rate. 
The resulting mean and standard deviation of logio ^ 
shown in figure 14, and we observe that there is essen- 
tially no correlation between colour ratio and turbulence. 
There is a hint that the standard deviation is slightly lower 
for the most extreme colour ratios; however this is because 
there are very few samples here. These statistics are con- 
sistent with the case study in section 4. Klett (1995) and 
Breon and Dubrulle (2004) performed theoretical calcu- 
lations of the influence on turbulence on crystal orienta- 
tion, and these also indicate that for the observed levels 
of £, crystals can be oriented to well within 1° of hori- 
zontal. Our observed range of e values are consistent with 
other ice cloud observations (eg. Gultepe and Starr 1995, 
Bouniol et al 2003). We note that the bulk of our e values 
are an order of magnitude lower than in Klett's calculation 
for 'weak' turbulence, and suggest that in most stratiform 
clouds, the ability of turbulent eddies to influence the air 
flow on the scales relevant to ice crystals is probably min- 
imal. 

The lack of correlation observed above leads us to 
believe that the colour ratio is primarily sensitive to the 
fraction of ice particles which are pristine planar crystals 
with Re close to 10. 

5.4 Supercooled layers 

As mentioned in section 3, it is our experience that super- 
cooled liquid layers are a frequent source of ice crystals 
producing specular reflection. Young et al. (2000), Piatt 
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Figure 14. Effect of turbulence on crystal orientation: line shows 
mean of logio ^ measured by the radar as a function of colour ratio; 
error bar shows 1 standard deviation in logio ^• 

(1977), Sassen (1984) also observed specular ice crystals 
falling beneath supercooled liquid layers. We have inves- 
tigated this link in a more quantitative manner by isolating 
liquid cloud using the off-vertical ceilometer in the man- 
ner described above, and then picking out the occasions 
when ice crysals are observed falling directly below this 
liquid layer. By the nature of the algorithm, we only pick 
out the lowest supercooled layer in multi-layered cloud 
systems. The fraction of these profiles affected by spec- 
ular reflection is plotted in figure 15 as a function of the 
temperature of the liquid layer. It is immediately appar- 
ent that specular reflection is very frequent in these cloud 
types, with 80% of ice falling below supercooled lay- 
ers between —12 and — 20°C affected. This is consistent 
with in-situ observations by Locatelli et al (1983), Hogan 
et al (2003a) and Field et al (2004) who all observed 
pristine vapour-grown planar crystals below supercooled 
layers in this temperature range; it is also consistent with 
our case study in section 4. Ryan et al (1976) found all 
crystals grown at these temperatures to be simple planar 
types. Such crystals therefore indeed seem to be charac- 
teristic of mid-level supercooled clouds as suggested by 
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Figure 15. Fraction of ice cloud observed below supercooled liquid 
which contains some specular reflection, binned as a function of 
the temperature of the liquid cloud. Horizontal lines indicate 95% 
confidence intervals for these fractions (assuming that if multiple 
samples were made from the same underlying population, the 
fractions derived from those samples are normally distributed). 



Hogan et al. and Field et al. - this may be a useful signa- 
ture for the possible presence of supercooled liquid water 
in optically thicker clouds where the liquid cannot be seen 
by the lidar, but the specular reflection from the ice crys- 
tals falling out is readily observed. 

There is a sharp decrease in the fraction of specular 
profiles at temperatures warmer than around -10°C; any 
crystals nucleated in such conditions are likely to be 
columns and needles, although there is still specular 
reflection observed in a significant fraction of cases. The 
lack of ice nuclei means that many of these warmer 
supercooled layers which are observed to be precipitating 
ice crystals are likely to be seeded from colder clouds 
above, in a similar way to the liquid layers embedded 
in the ice virga in fig 2. At temperatures colder than 
^ — 35°C there are no cases in our 17 month data set, 
which seems sensible as this is close to the limit for 
homogeneous freezing. 

Even supercooled layers at cold temperatures < 
— 22°C often have specular reflection in the crystals falling 
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Figure 16. Filled contour plot showing distribution of colour ratios 
in specularly reflecting ice cloud falling beneath supercooled liquid 
layers. Grey scale is logio (frequency) in the 2.5° C x 0.5dB bins. 
Contours range from 1.0 to 2.5 in steps of 0.3. 



below, around two- thirds of the time. This conclusion 
is also supported by Heymsfield et aUs (1991) observa- 
tion of a highly supercooled altocumulus layer — 30°C) 
with specular ice crystals precipitating out. To investigate 
this behaviour in more detail, figure 16 shows the distri- 
bution of colour ratios measured in the specular ice, as a 
function of supercooled layer temperature. This distribu- 
tion demonstrates that specular reflection is strongest in 
the ice falling from supercooled layers which lie between 
— 10 and — 20°C (where the droplets will freeze into planar 
monocrystals eg. Ryan et al. 1976). At temperatures below 
— 20°C, the magnitude of the colour ratio decreases sub- 
stantially, indicating that planar monocrystals are typically 
a minority in these colder clouds. In what follows, we 
speculate on the origin of the specular crystals which are 
present in cold clouds and why they are fewer in number 
than at warmer temperatures, by linking previous labora- 
tory and observational measurements of crystal growth to 
the data in figures 15 and 16. 

Growth at temperatures < — 22°C has traditionally 
been associated with columnar crystal habits (Pruppacher 
and Klett 1997), and one may hypothesise that the 



QJRMS 



QJRMS 00: 1-26 (2009) 



20 C. D. WESTI 

observed specular reflection is the result of Parry-oriented 
columns (major axis horizontal, prism face horizontal), 
rather than plates; such a case has been documented by 
Sassen and Takano (2000). However observations of the 
associated Parry arc are very rare, and Sassen and Benson 
(2001) have argued that in general only planar ice crystals 
will produce specular reflection. 

Perhaps a more likely explaination of our observa- 
tions is that planar crystal growth is occurring at colder 
temperatures than expected. Planar crystals have been 
observed to grow at cold temperatures in a number of 
studies. Aufm Kampe et al. (1951) grew mainly plate- 
like polycrystals from a supercooled water cloud in a cold 
chamber at — 22°C and — 29°C, along with a smaller frac- 
tion of simple plates and columns. Bailey and Hallett 
(2002,2004) nucleated and grew ice crystals on kaolin- 
ite and bare glass fibres, and found that plate-like growth 
in mono- and poly-crystalline forms dominated down 
to approximately — 38°C. Bacon et al. (2003) nucle- 
ated and observed the very early stages of growth of 
supercooled droplets suspended in an electrodynamic trap 
and observed plate-like growth at temperatures as cold 
as — 38°C. Fleishauer et al. (2002) observed a mixed- 
phase cloud topped at — 30°C: the ice particle images 
from that case study seem to show a mixture of pla- 
nar mono- and poly-crystals. Sakurai (1968) observed ice 
crystals in supercooled fogs between —20 and -28°C and 
found irregular polycrystals, thick plates and short (almost 
1:1 aspect ratio) columns. These laboratory and observa- 
tional data suggest that crystals nucleated and grown from 
droplets at these cold temperatures often favour plate-like 
growth, albeit typically as a poly crystal rather than simple 
monocrystalline plates. 

Pitter and Pruppacher (1973) have investigated the 
freezing of fairly large liquid water droplets (180-800/im 
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in diameter) by immersion and contact nucleation (using 
clay particles), and observed that for a given supercooling 
AT, there was a critical size for the droplet to grow into 
a poly crystal: ex (AT)"^. Bacon et al. (2003) found 
that their 40/im diameter droplets grew into polycrystals 
at temperatures colder than ^ -22°C. Based on this, one 
might expect a smaller 20/im droplet to develop into a 
polycrystal at ^ -27°C; for a 10/im droplet, monocrystals 
would dominate down to ^ — 35°C. These thresholds are 
approximately consistent with the extrapolation of Pitter 
and Pruppacher' s data to cloud droplet sizes. 

We suggest therefore that our observed specular 
reflection in cold clouds is the result of small cloud 
droplets freezing into simple (thick) plates which pro- 
duce specular reflection, whilst the larger droplets in 
the droplet distribution freeze into irregular polycrystals. 
Whereas in clouds warmer than — 20°C the vast majority 
of cloud droplets will freeze as monocrystals irrespective 
of droplet size, for colder clouds an increasing fraction of 
the droplets frozen will be polycrystalline, and incapable 
of producing specular reflection. 

We note that although on average specular reflection 
in cold clouds is much weaker, there are still individual 
cases where it is large -lOdB) - these may be exam- 
ples where the droplets are particularly small. For the 
low liquid water content stratiform clouds which dominate 
our observations, we anticipate there may be a stronger 
weighting towards smaller droplets (and therefore more 
single crystals) than in more convective-type clouds. Note 
that the droplets in Heymsfield et al. 's (1991) altocumu- 
lus layer with specular reflection in the virga below had a 
mean diameter of only ^ 10/im. 

5.5 Layers of planar crystals 

Our observations (eg. section 3) indicate that specular 
crystals form widespread layers stretching over many 
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kilometres; this also fits in with previous fidar obser- 
vations (Piatt 1977). Aircraft observations by Korolev 
et al. (2000) suggest that branched planar crystals often 
occur in 'cells' many kilometres across. To get an idea of 
the typical scale of these planar crystal layers, we have 
identified areas of cloud which have at least one common 
range gate with a negative colour ratio in two consecu- 
tive vertical profiles - in this way we are able to identify 
continuous layers of specular reflection, and characterise 
their horizontal extent. Of course there is the complica- 
tion that the lidar is frequently blocked by low level liquid 
clouds. To ameliorate this effect (which will tend to bias 
our estimates too low) we allow gaps of up to 5 minutes 
in the layer so long as there is still an adjacent region of 
specular reflection on both sides of the gap. Using the 
wind speeds from the Met Office forecast model, figure 
17 shows the mean horizontal extent of the planar crys- 
tal layers as a function of the average temperature of the 
layer. It is apparent that in mid-level clouds the average 
extent is approximately 15 km, whilst at colder temper- 
atures the layers appears to be less long-lived with an 
average extent of approximately half that value. In light of 
the discussion above regarding specular reflection below 
liquid layers, this may perhaps be the result of the more 
transient nature of supercooled liquid clouds at these cold 
temperatures, which typically are only ~ 5 km in width at 
— 35°C (Hogan et al. 2003b). There is a also a systematic 
decrease in layer extent at warm temperatures in figure 
17. The reason for this is unclear, although it may repre- 
sent a shift in weighting of the statistics from widespread 
supercooled altocumulus virga to the base of deeper ice 
clouds associated with approaching fronts. We note that 
at all temperatures the distribution of layer widths is quite 
broad, with some cases lasting for hours on end, whilst 
others only occur in isolated pockets or fall streaks, and 




5 10 15 20 25 

Average horizontal extent of specular layers [km] 

Figure 17. Mean horizontal extent of specularly reflecting layers, 
as a function of the average layer temperature. Error bars indicate 
uncertainty in the mean values. 



we observed layers from ~ 1 km up to more than 100 km 
in extent. This large variability is also consistent with 
Korolev et a/.'s observations. 

The depth of the specular layers has also been char- 
acterised - we estimate a mean depth of 600m, approxi- 
mately independent of temperature. This figure is a con- 
volution of the true layer depth and the extinction of the 
lidar beam in optically thick clouds, and is likely an under- 
estimate. 

5.6 B ackscatter and radar reflectivity 

Figure 18 shows distributions of attenuated backscatter 
recorded from each lidar, with separate bars indicating 
specular and normal ice cloud. The data from the verti- 
cally pointing 1.5/im lidar shows that the backscatter from 
clouds affected by specular reflection is typically a factor 
of 5dB larger than that recorded for normal clouds. From 
the 905nm data we observe the opposite, ie. clouds where 
specular reflection is present typically have a backscatter 
which is typically around 5dB lower than that recorded in 
normal ice cloud. This suggests that cloud where specu- 
lar reflection is occurring may tend to be optically thinner 
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Figure 18. PDFs of lidar backscatter recorded by the vertical 
1.5 fim and off- vertical 905 nm lidars. Black bars are the PDFs for 
cloud identified as containing specularly reflecting ice crystals, grey 
bars are PDFs for normal ice cloud. 



than average. However, the interpretation of this result is 
subtle since the lidar is attenuated as it passes through the 
cloud, and our experience suggests that specular reflection 
is rarely observed at cloud base (because the crystal facets 
become rounded as they evaporate), which may bias the 
statistics towards lower values of off-zenith backscatter. 

Radar reflectivities have been measured from the 
3 5 GHz cloud radar at Chilbolton, and data were analysed 
for the 7 month period June 2007 to January 2008. Figure 
19 shows the distribution of radar reflectivity measured 
in ice clouds which are also observed by both lidars, 
and it appears that the reflectivity distribution in clouds 
which are affected by specular reflection (black bars) is 
almost identical to the reflectivity distribution for all ice 
cloud (grey lines). The suggestion then is that the ice 
water content in cloud where there is specular reflection 
is representative of the complete subset of ice cloud 
observed by both lidars. 
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Figure 19. PDFs of radar reflectivity for clouds containing specu- 
larly reflecting crystals (black bars), all ice cloud detected by both 
lidars (grey lines). 

6 Summary 

'Mirror- like' specular reflections from horizontally ori- 
ented planar ice crystals have been observed by compar- 
ing measurements from a vertically pointing Doppler lidar 
and a second lidar pointing slightly off-zenith. This com- 
parison allows supercooled liquid cloud (which is highly 
reflective at both elevations) to be separated from layers of 
oriented planar crystals (which only produce a strong sig- 
nature at vertical). Whilst most previous studies of spec- 
ular reflection have focus sed on glaciated cloud profiles, 
we find that strong specular reflection is most common 
in mid-level mixed-phase clouds, and is observed 85% of 
the time in ice falling from supercooled layers at -15°C. 
These crystals typically formed well defined layers with 
a mean extent of 15km, although numerous cases 100km 
or more across have been observed. Similar vast regions 
of planar ice crystals have been observed beneath super- 
cooled boundary-layer clouds in Alaska (A. Korolev, per- 
sonal communication). 

A case study of a long-lived supercooled layer with 
planar ice crystals precipitating below was presented as a 
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Future work will focus on the analysis of data from a 
depolarisation channel added to the Doppler lidar in 2008, 
providing high resolution measurements of depolarised 
backscatter and depolarised Doppler velocity. 
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typical scenario where strong specular reflection is likely measurements of light pillars by Sassen (1980). The asso- 

to be observed, and the backscatter at vertical was a factor ciated crystal dimensions are approximately 400-1200/im 

of 30 (15dB) larger than observed off-zenith. Simultane- depending on the assumed crystal habit. 

ous polarimetric radar measurements confirmed that the 

general crystal population was well oriented with values 

of Zdr typically ^ 2dB. Observations at the base of the 

ice virga reveal that specular reflection does not occur in 

this region, even though the Zdr measurements show that 

the crystals are still well oriented, and this is evidence that t . i • . i £i 

Interestingly, significant specular reflection is 

the crystals are evaporating rapidly and becoming rounded , a ■ € ^^^^ ^ iji -ji 

& ^ o observed m ice railing below supercooled liquid layers at 

in shape. This may be a useful observational technique , , oojoro-oi- • ^^ 

^ ^ temperatures between —22 and — 35 C. Freezing or large 

to identify whether crystals are evaporating layer-by-layer ijiw. . .-n j 

& ^ ^ ^ cloud droplets at such temperatures typically produces 

at small subsaturations, or whether the evaporation rate is . , i.ii i .11 1 .1 

^ irregular plate-like poly crystals; however poly crystal 

rapid and uncontrolled in drier air (Nelson 1998, Hallett .1 , x-.- r ^ ^ ^ ■ ^ 

^ development is a strong function or droplet size, and 

et al. 2001). suggest that our observations can be explained by 

freezing of relatively small droplets, which develop into 

simple thick plates and produce specular reflection. This 

may be more favourable in the gently forced stratiform 

clouds that we are sampling; in more convective clouds 

with bigger droplets polycrystals would almost certainly 

dominate. 



High temporal resolution Doppler radar measure- 
ments reveal that there was appreciable convective over- 
turning at cloud top due to radiative cooling of the liquid 
layer to space, and this turbulence extended down into 
the ice virga. Despite this turbulent air motion (e 3 x 
10~^m^s~^), strong specular reflection was still observed 

in this region of virga, indicating that turbulence does Finally, whilst we have observed that specular reflec- 
not have a strong effect on crystal orientation in strati- tion is strong in many cases, we see from figure 1 1 that in 
form clouds. This is backed up by the lack of correlation the majority of cloud it occurs at a relatively weak level 
between the measured colour ratio and e in section 5.3. with colour ratios close to dB, and indeed at all temper- 
atures the modal colour ratio is positive (ie. no specular 
Lidar Doppler measurements indicate that these hor- reflection). This lends support to the in-situ observations 
izontally oriented crystals fall slowly, with a mean veloc- of Korolev et al (2000) which indicate that most ice par- 
ity of ^ 0.3ms~^ and a weak temperature dependence, in tides in thick stratiform clouds are 'irregular', probably a 
contrast to normal ice cloud where the observed sedimen- mixture of complex polycrystals and aggregated particles 
tation velocities in warm clouds are twice as fast as those (Stoelinga et al 2007), but with few planar monocrystals 
at -40° C. Using velocity-diameter relationships, we esti- mixed in. The important exception to this rule appears to 
mate that the most strongly oriented planar crystals have be thin supercooled layer clouds, particularly those in the 
a Reynolds number in the region of Re ~ 10 with a fac- range —10 to — 20°C, which our measurements indicate 
tor of 2 uncertainty. This value is consistent with previous are a prolific source of pristine planar monocrystals. 
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Appendix A: Calibration of the 1.5/im Doppler lidar 

We follow O'Connor et al. (2004) and compute the extinc- 
tion to backscatter ratio S for a gamma distribution of 
liquid water droplets for the new lidar wavelength of 
1.5/im. This is plotted as a function of the median volume 
droplet diameter I^o for various values of the dispersion 
parameter ji in figure 20, where the droplet diameter dis- 
tribution takes the form n{D) ex (i:>/i:>o)^ exp[-(3.68 + 
IJi)D/ Dq] . For realistic values of the median volume diam- 
eter Do (8-40/im) and /i (2-10), 5 = 20 ± 20%. For an 
optically thick liquid cloud which completely attenuates 
the lidar signal the integrated backscatter through the 
cloud is simply {2S)~^ (see O'Connor et al.) and the cal- 
ibration is applied simply by scaling the backscatter so 
that the integrated backscatter measured in thick stratocu- 
mulus cloud matches S = 20. Because of the very narrow 
field of view, multiple scattering is neglected. 

Appendix B: Differential reflectivity calculations 

We estimate Z^r at 45° elevation for horizontal planar 
crystals by approximating them as oblate spheroids com- 
posed of a homogeneous mixture of air and ice inclusions. 
The incident wave is alternately transmitted in two orthog- 
onal polarisations: one parallel to the horizontal plane 
('H'), and a second at 45° from the horizontal plane ('V'). 
Both are orthogonal to the direction of propagation and 



Extinction/Backscatter [sr] for |Li=2-10 




Dq [|Lim] 



Figure 20. Extinction to backscatter ratio for a gamma distribution 
of liquid droplets with /x = 2 (darkest gray) to 10 (lightest gray). 



have amplitude Eq. Resolving the incident wave along the 
spheroid's major and minor axes the amplitudes are: 

^0 - along major axis 1 (parallel to H) (2) 
E^l — along major axis 2 (± to H) (3) 
Eq/^J2 - along minor axis (4) 

inducing dipole moments: 



or 

0.5^0 (<^maj or + <^minor) 



in the H channel (5) 
in the V channel (6) 



where 



:major, minor 



are the polarisabilities: 



ai = particle volume x 



e- 1 



1 + L,(e-1)* 



(V) 



The geometric factors Li for spheroids are detailed in 
Bohren and Huffman (1983); the permittivity of the air 
ice mixture is calculated using Maxwell-Garnett's rule 
for spherical inclusions (see Bohren and Huffman) - this 
prescription depends only on the volume fraction of ice in 
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the mixture fy = (volume of ice) /(volume of spheroid). 
The differential reflectivity is then simply: 



Zdr = 10 logio 



a 



major 



(8) 



which is maximised for thin plates of solid ice (/^ = 1) at 
Zdr = 3.6dB. Hogan et al. (2002) have performed more 
extensive calculations of this parameter for 0° elevation. 
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